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A delicate interplay between the anomalous magnetic moments of the proton 
and neutron makes, in magnetic fields B > 2 x 10^^ T, the neutron stable and for 
fields B > 5 X 10^^ T the proton becomes unstable to a decay into a neutron via /3 
emission. Limits on the field strengths for which these arguments hold are presented 
and are related to questions of vacuum stability in the presence of such fields. Possible 
astrophysical consequences are discussed. 
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I. INTRODUCTION 



Very intense magnetic fields have been postulated to exist in connection with some as- 
trophysical objects. Fields with strengths larger than 10^^ T are associated with supercon- 
ducting cosmic strings and recently a proposal to explain extragalactic gamma ray bursts 
involves fields of around 10^'^ T [Q. We have been investigating several questions related to 
the existence of such fields and to the behavior of elementary and composite states in such 
environments. In a previous paper we discussed the possible breakdown of constant 
magnetic fields with strengths beyond 10^^ T. Other mechanisms that destabilize strong 
magnetic fields have been proposed @], and here we add another intriguing phenomenon to 
this list. For fields greater than 5 x 10^^ T the proton becomes heavier than the neutron 
and decays into the latter by positron emission. 

In Section II we study the behavior of an proton, neutron and electron in an intense 
magnetic field and find the aforementioned amusing result that the proton becomes unstable 
against neutron, positron and neutrino decay. The decay rates and spectrum are obtained 
in Section III. Whereas the discussion of the behavior of the electron is on firm footing, 
questions may be raised as to the validity of our treatment of the proton and neutron; these 
questions and the stability of the vacuum in the presence of strongly interacting particles 
with anomalous magnetic moments of non-electromagnetic origin are discussed in Section 
IV. Conclusions and experimental consequences are presented in the last section. 

II. LOW LYING STATES FOR PARTICLES IN UNIFORM MAGNETIC FIELDS 

The quantum mechanics of a Dirac particle with no anomalous magnetic moment in a 
uniform external magnetic field is straightforward. We shall present the results for the case 
where particles do have such anomalous moments. In reality, in fields so strong that the 
mass shifts induced by such fields are of the order of the mass itself one cannot define a 
magnetic moment as the energies are no longer linear in the external field. Schwinger ^ 
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calculated the self energy of an electron in an external field and we shall use his results 
subsequently. We cannot follow this procedure for the proton or neutron as we do not have 
a good field theory calculation of the magnetic moments of these particles, even for small 
magnetic fields; all we have at hand is a phenomenological anomalous magnetic moment. 
However, for fields that change the energies of these particles by only a few percent, we will 
consider these as point particle with the given anomalous moments. In the Section IV we 
will discuss possible limitations of this approach. 



A. Protons in an External Field 



The Dirac Hamiltonian for a proton with a uniform external magnetic field B is 



H = Qi-{p- eA(r)) + pMp 



^ - 1 ) /3S ■ B . 



(1) 



The vector potential A(r) is related to the field by A(r) = |r x B and gp = 5.58 is the 
proton's Lande g factor. We first solve this equation for the case where the momentum 
along the magnetic field direction is zero and then boost along that direction till we obtain 
the desired momentum. For B along the z direction and Pz = the energy levels are 



En' 



2eB{n + ^) - e5s + Mp^ 



9p 



1 fis, 



(2) 



In the above, n denotes the Landau level, m the orbital angular momentum about the 
magnetic field direction and s = ±1 indicates whether the spin is along or opposed to that 
direction; the levels are degenerate in m. n = and s = +1 yield the lowest energy 



E = M„ = M„- 



IB. 



(3) 



As we shall be interested in these states only we will drop the n and s quantum numbers. 
The Dirac wave function for this state is 



^m,p,=o(r) 








(fornix, y) 



(4) 



are the standard wave functions of the lowest Landau level; 



(t>m{x-,y) = 



m+l 



\/wrn}. 



-[x + iy]"*exp 



Boosting to a finite value of Pz is straightforward; we obtain 



with a wave function 



^ cosh 6* ^ 





sinh^ 




ol-PzZ 



where 29, the rapidity, is obtained from tanh2^ = p^/ Em{pz)- 
In the non-relativistic limit the energy becomes 
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Pz 



2M 



and the wave function reduces to 



27r' 



(5) 



(6) 



(7) 



(8) 



(9) 



B. Neutrons in an External Field 



For a neutron the Dirac Hamiltonian is somewhat simpler 



e fg 



- ) /3S • B . 



(10) 



2M„ V 2 y 

with gn = —3.82. Again for Pz = the states of lowest energy, the ones we shall be interested 
in, have energies 
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^(P±,P. = 0) 



2M„ V 2 



(11) 



Boosting to a finite Pz we obtain 



(12) 



The wave functions corresponding to tfiis energy are 



= -—^u{p,s = -1) 

[ZTT) 2 



(13) 



wfiere u{p,s = —1) is tlie standard spinor for a particle with momentum p, energy 



p2 + (not the energy of Eq. ([T2|)) and spin down. 
In the non-relativistic hmit 



2M„, 



(14) 



and the wave functions are 



^p(r) 



1 




(27r)§ 



(15) 



C. Electrons in an External Field 

We might be tempted to use, for the electron, the formalism used for the proton with 
the Lande factor replaced hj Qe = 2 + a/vr. However as we shall see for magnetic fields 
sufficiently strong as to make the proton heavier than the neutron the change in energy 
of the electron would appear to be larger than the mass of the electron itself. The point 
particle formalism breaks down and we have solve QED, to one loop, in a strong magnetic 
field; fortunately this problem was treated by Schwinger ||^. The energy of an electron with 
Pz = 0, spin up and in the lowest Landau level is 
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Em,p^=0 — Me 



a , f2eB\ 
IH In ' ' 



(16) 



27T \Me^ I _ 

For field strengths of subsequent interest tliis correction is negligible; the energy of an 
electron in the lowest Landau level, with spin down and a momentum of is 



+ M, 



(17) 



and with wave function similar to those of the proton 



V'm,p,(r) 



' ^ 



cosh 9 




sinh 6 



(18) 



where the boost rapidity, 26*, is defined bellow Eq. (^ while the Landau level wave function 
is defined in Eq. (^. The reason the complex conjugate wave function appears is that the 
electron charge is opposite to that of the proton. 

D. Decay Kinematics 



From Eq. (^) and Eq. (p!^ ) we note that the neutron becomes stable against /5-decay 
when the following inequality is satisfied 



2M„ V 2 



2M„ V 2 ; - « p 



(19) 



or i? > 2 X 10^^ T. On the other hand the proton becomes unstable for decay into a neutron 
and a positron whenever 



2M„ V 2 / 2Mp \ 2 



^ - 1 ) 5 ~ 0.12^^5 > M„ + Me -Mp, 



(20) 



or for 5 > 5 X 10^^ T. We shall now turn to a calculation of the life time of the proton in 
fields satisfying this inequality. 
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III. PROTON LIFE TIME 



A. Proton, Neutron and Electron Fields 



With the wave functions of the various particles in the magnetic fields we may define field 
operators for these particles. For the proton and electron we shall restrict the summation 
over states to the lowest Landau levels with spin up, down respectively; for magnetic fields 
of interest the other states will not contribute to the calculation of decay properties. For the 
same reason, the neutron field will be restricted to spin down only. The proton and neutron 
kinematics will be taken as non-relativistic. 



^0^ 





1 



-ipzZ 



(21) 



with 4>m{x, y) defined in Eq. (||) and the energy, Emipz) in Eq- (i)- CLm{Pz) is the annihilation 
operator for a proton with momentum p^z and angular momentum m; bm{Pz) is the same 
for the negative energy states. For the neutron the field is 



a{p) 



1 




oipr 



{271 








-ipr 



{2n 



(22) 



with an obvious definition of the annihilation operators. For the electron we use fully 
relativistic kinematics and the field is 
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Mr) ^ 



amiPz) 



' ^ 

cosh Q 


sinh Q 

\ 





cosh 6' 


sinh 6 



-iPzZ 



(23) 



B. Decay Rates and Spectrum 

The part of the weak Hamiltonian responsible for the decay p ^ n + + Ue is 

H = ^J d'x^^j^il + 75)vl/,'l',7'^(l + 75)v|/, . (24) 

For non-relativistic heavy particles the matrix element of this Hamiltonian between a proton 
with quantum numbers = , m = m^, a neutron with momentum a neutrino with 
momentum pj, and an electron in state m = rrif and with Pz^e is 



(H) = ^ Me./2)^—6{pz,+Pz,.+Pz,^) 

J dx dy (j)*^^ (x, y)(j)m, (x, y) exp[-i(px,ri + P±,u) ■ r^] ; (25) 

6i, is the azimuthal angle of the neutrino. The integral in the above expression can be 
evaluated in a multipole expansion. Note that the natural extent of the integral in the 
transverse direction is 1/^/eB whereas the neutron momenta are, from Eq. (^), of the 



order of y/0.12eB] thus setting the exponential term in this integral equal to one will yield 
a good estimate for the rate and spectrum of this decay. The positron spectrum is given by 



rfr _ 4 Gj^Mp Ee + Pz 



(A - EeY ; (26) 



dpz^e 3 (27r)6 E, 
where A = 0.12/ijv-B — M„ + Mp. For A ^ Mg the total rate is easily obtained 
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r 



3 (27r)6 



(27) 



The lifetime is r ~ 1.5 x 1(}^{W^^T/By s. 

IV. FIELD STABILITY IN THE PRESENCE OF COMPOSITE HADRONS 

In Refs. it was noted that vector particles with anomalous magnetic moments induce 
instabilities for large magnetic fields. Such a mechanism is not induced by spin one-half 
fields, even with anomalous moments. It is however clear that the discussion of the previous 
section regarding protons and neutrons has to break down at a sufficiently strong field. Let 
us look specifically at the neutron case. From Eq. (|12D we see that for M„ — fiB < positive 
and negative energy levels cross and pair creation becomes possible. Following the analysis 
of Ref. we can show that the vacuum does not decay into neutron-antineutron pair. (We 
cannot rule out the possibility that the strong interactions responsible for the anomalous 
moment could induce, for example, pion pair creation.) What occurs, is that the vacuum 
becomes a linear combination of the old vacuum and two fermion-antifermion pairs; a mixing 
with a single fermion-antifermion pair is ruled out on parity grounds. 

It is amusing to look at this problem in two dimensions^. A Lagrangian for a charged 
Fermion with an anomalous moment is 



In two dimensions there is only one field component, F^'^ = E. We note that for non-zero 
charge, e, a true instability develops, whereas for e = and non-zero fi the field (p acquires 
a finite expectation value, corresponding to fermion pairs in the original model. 



^This two dimensional argument is due to A. Schwimmer. 



(28) 



The bosonized version of this Lagrangian is 




(29) 
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From the above we expect the arguments of the previous sections to break down for 
critical fields, ~ M^/e ~ 10^^ T. However, the composite nature of the nucleons will 
cause a breakdown for lower fields. A reasonable estimate would be a critical field for the 
constituent quarks, B ~ M^/e ~ 10^^ T. Thus our analysis should be valid for fields below 
this value and that leaves a window, 5 x 10^^ T < 5 < 10^^ T where we expect the proton 
to be heavier than the neutron and to decay into it. 

V. CONCLUSIONS AND EXPERIMENTAL CONSEQUENCES 

We studied the mass evolution of protons, neutrons and electrons in strong magnetic 
field and concluded that the proton will decay /5 decay, p ^ n + e'^ + Ue, in a sufficiently 
strong external magnetic field. The new mechanism for positron generation might have 
astrophysical interest. There is indeed an overabundance of positrons as compared with 
that accounted for by existing mechanisms [Q. 
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